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Abstract 

We investigate the effect that Rashba spin-orbit coupling has on 
the low energy behaviour of a two dimensional magnetic impurity 
system. It is shown that the Kondo effect, the screening of the mag- 
netic impurity at temperatures T < Tk, is robust against such spin- 
orbit coupling, despite the fact that the spin of the conduction elec- 
trons is no longer a conserved quantity. A proposal is made for how 
the spin-orbit coupling may change the value of the Kondo tempera- 
ture Tk in such systems and the prospects of measuring this change 
are discussed. We conclude that many of the assumptions made in 
our analysis invalidate our results as applied to recent experiments 
in semi-conductor quantum dots but may apply to measurements 
made with magnetic atoms placed on metallic surfaces. 

Key Words: Kondo effect; Rashba spin-orbit coupling; mag- 
netic impurities; strongly-correlated electrons 

1 Introduction 

The physical behaviour of a single magnetic impurity interacting with a 
large number of conduction electrons has been, and continues to be, a 
fascinating subject of inquiry for over four decades. The emergence of 
a strong correlation between the impurity and the conduction electrons 
at low temperatures, first predicted by Jun Kondo [TJ, gives rise to the 
well-known Kondo effect wherein the impurity is effectively screened by 
forming a singlet with the surrounding electrons. Much still remains to 
be understood regarding strongly correlated quantum many-body physics. 
The Kondo effect, as a well understood example of a strongly correlated 
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phenomenon, has proven to be a fertile topic for further study in this 
important field [2]. 

In recent years, the study of these magnetic impurity systems has seen a 
resurgence due to the great advances in experimental technology that allow 
one to construct, control, and manipulate objects on the nanometer scale. 
In particular, the Kondo effect has been observed in numerous experiments 
involving quantum dots constructed in semiconductor heterostructures [3] , 
as well as in systems where a magnetic atom is placed on the surface of 
metal g]. 

Both of these systems involve a local moment, either that of an odd 
number of localized electrons (as in the quantum dot) or of a partially filled 
outer shell of an atom interacting with mobile electrons confined to move 
in two dimensions. The confining potential breaks the inversion symmetry 
that is manifest in the conduction electrons and gives rise to a coupling 
between their spin and momentum [5]. Such spin-orbit coupling (coined 
as Rashba spin-orbit coupling to differentiate it from coupling arising from 
other mechanisms which may break the inversion symmetry) has been an- 
alyzed theoretically and observed experimentally in both semiconductor 
quantum wells [7] and in metallic surface states [5J OH [TO] ■ 

The presence of such a spin-orbit coupling means that the SU (2) sym- 
metry of the conductions electrons is no longer present and that the spin 
of the electrons is no longer an eigenstate of the system. Given this fact, it 
seems, at first, remarkable that the Kondo effect, which depends crucially 
on the magnetic exchange interaction between the local moment and the 
spin of the conduction electrons, is observed at all in such two dimensional 
systems. This issue was first examined in where it was shown that 
the presence of time reversal symmetry in the Kondo interaction term pre- 
serves the spin coherence of the electron propagator. That is, although 
spin is no longer a "good" quantum number, the scattering trajectories of 
the electron that would destroy its spin coherence in multiple interactions 
with the impurity are exactly cancelled by the time reversed paths. While 
such a powerful argument may explain why the Kondo effect survives the 
apparent breaking of spin coherence, it is not precise enough to determine 
the change in the low energy behaviour of such impurity systems in the 
presence of a spin-orbit coupling in general. A more detailed calculation is 
required in order to gain such information. 

There have been several studies that have investigated the role of spin- 
orbit effects in specific semiconductor mesoscopic devices where the Kondo 
effect is manifest [HI [13l [TH [15] . In these studies, the spin-orbit interac- 
tion either takes place within a quantum dot coupled to conducting leads 
or within one dimensional wires, often coupled to a localized magnetic mo- 
ment. In this paper, we derive a model to study the effect that the Rashba 
spin-orbit coupling has on the low energy behaviour of a two dimensional 
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system with a magnetic impurity. The generic nature of the derived model 
may not be as readily applicable to the analysis of a specific experiment 
as those cited above, but it does have the advantage of elucidating some 
general features of the physical interplay between the Kondo and Rashba 
interactions that might be less transparent in more complicated quantum 
dot studies. 

The model is presented in S}2]and an effective model relevent to the study 
of the low energy behaviour is derived. The analysis of the effective model 
suggests that a more physical property of the system should be calculated 
in order to get a better sense of how the spin-orbit coupling influences 
its behaviour. In S£2 the resistivity of the full two dimensional model is 
calculated and compared to Kondo's classic result pQ. The combination of 
this result together with an analysis of the effective model leads to various 
conclusions and predictions which are discussed in 21 

Throughout the analysis we use the following conventions: vector quan- 
tities are represented by a boldface character and the same character with- 
out boldface denotes the vector's magnitude. Repeated Greek indices are 
assumed to be summed over unless one of them is enclosed in brackets. 
Finally, we work in units in which h = 1 except when estimating the values 
of certain quantities in which case the precise units will be given. 

2 Model Hamiltonians 

2.1 Full Two Dimensional Model 

We first consider electrons confined to travel in the xy plane as the result 
of some non-constant electric potential in the z direction. The motion of 
the electrons in the presence of this electric field gives rise to a coupling 
between the spin of the electron and its momentum, which can be expressed 
by the Rashba Hamiltonian [5, 6 J 

Our notation is as follows: m is the effective mass of the electron, 'J/j^ 
is the operator that annihilates an electron of momentum k = (k x ,k v ) 
and spin fi and satisfies the anticommutation relationship {'S'jL, 'i'k'u} — 
5^5^ (ft, ft'), a is the strength of the Rashba coupling which is propor- 
tional to the derivative of the electric potential in the z direction, and cr* 

are the Pauli matrices so that \I '! ir^f is twice the i th component of 
the spin of the electron. 
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This Hamiltonian can easily be diagonalized to be of the form 

Hli = T^yl d2k J24Al a A ka . (2) 
The Rashba bands are non-degenerate and take the form 

4 = ^(k±k R ) 2 -e R (3) 

where we have defined the Rashba momentum scale k R :— ma and the cor- 
responding Rashba energy scale er := k R /2m. The relationship between 
the operators that create particles in these Rashba bands and the original 
electron operators is 

where 9 is the angle of the momentum k in polar coordinates and the 
normalization has been chosen such that {4„, ^Wb} = SabS^ {k, k'). 

Now consider a spin-1/2 magnetic impuritjo at the origin, described by 
the spin operator S. We write the interaction between this spin and the 
electrons using the so-called s-d exchange Hamiltonian or Kondo Hamilto- 
nian 

Hk = (2^ J d2kd2k ' Jkk ' s ' ^L^f*^' ( 5 ) 

where V is the two dimensional volume of the system. This interaction 
can be thought of as arising from an Anderson impurity model [17] of a 
two dimensional Fermi gas interacting with a localized electronic state with 
strong Coulomb repulsion. The details of this relationship (in the absence 
of spin-orbit coupling) are discussed in appendix [Al where the specific form 
for Jfcfc' is derived for various models of the hybridization between the local 
state and the electronic Fermi gas. 

In the absence of Rashba spin-orbit coupling, the physics of the Kondo 
Hamiltonian is very well known [5] . In the case of a ferromagnetic coupling 
( Jkk' < 0) between the impurity and the electrons, the impurity essentially 
decouples from the electrons at low temperatures. If the exchange interac- 
tion is antiferromagnetic, then there exists a temperature scale Tk, the so- 
called Kondo temperature, below which the magnetic impurity is screened 
by forming a singlet with the surrounding electrons. Many low-energy 
properties can then be derived (using a variety of techniques depending on 
the temperature regime) and most can be shown to depend on universal 

It is straightforward to generalize the results of this paper for the case of an impurity 
with higher spin. 



4 



scaling functions of T/Tk [2]- For the case where the interaction occurs at 
a single point, Jkk' = J and the Kondo temperature can be expressed to 
lowest order in J as 

T K = De~^ (6) 

where po is the density of states evaluated at the Fermi energy and D is 
an effective bandwidth. 

In semiconductor quantum dots, the value of the Kondo temperature is 
tunable though the maximum value is set by the size of the dot: the smaller 
the dot, the larger the temperature. In this way, Kondo temperatures as 
high as 1 K have been obtained pJ] . The Kondo temperature of a magnetic 
atom on a metallic surface can span quite a wide range depending on the 
nature of the impurity atom and which metallic surface is used. In the case 
of a Co atom on Au(lll), a Kondo temperature of Tjc ~ 70 K has been 
measured [ID) . 

To describe a system with both a Kondo impurity and Rashba spin-orbit 
coupling, we can use eq. (0| to write the Kondo interaction compactly in 
the Rashba operator basis as 

Hl< = 1(2^ / rf2fcd2fc ' Jkk ' S S kfe'4«A'& (7) 

where 

S& := iae ie 'S- - ibe~ ie S + + (l - abe- l(e - e '^ S z (8) 

and S := S x ±iS v are the raising and lowering operators for the impurity 
spin. On the right hand side of this equation, a, b = ±1 when the labels 
a, b = ± respectively. 

The full two dimensional Hamiltonian, which we call the Kondo- Rashba 
model herein, is then H = H R + H K where H R is given by eq. ^ and H K is 
given by eq. l[7]). Given the complicated nature of the interaction between 
the Rashba quasi-particles with non-degenerate bands and the magnetic 
impurity, it is not clear at this stage whether or not the well-known physics 
of the Kondo effect, briefly described above, is still manifest in the Kondo- 
Rashba model. 

2.2 Low Energy Effective One Dimensional Model 

The simplest form for the exchange interaction is to consider that it occurs 
at a single point. As discussed in appendix [A) this means that Jkk' = 
J. In this case, decomposing the operators into their different angular 
components reveals that only a finite number of angular modes couple to 
the impurity. 
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We first investigate the Kondo interaction in the absence of spin-oribt 
coupling. With k and 9 being the amplitude and angle of the momentum 
k, we can write 



* feM = ^=L= V ip kfim e imd (9) 
and using the orthonormality of (l/v / 27r)e l ™ e this implies that 

VVm = \/f 1^ dd e-^ 6 ^. (10) 



We have included the prefactors of Vk here so that the polar components tp 
of W behave as ordinary Fermionic operators {i> k ^ ml ipk'vn} — S^ l/ 5 mn S(k — 
k'). In terms of these angular components, the Kondo interaction ([5]) can 
be written as 

Hk = -r fn ,o / dkdk'Vkk 1 S-il)l Q tT llv il) k , V Q. (11) 
2 ( 27r ) Jo 

That is, the Kondo interaction is effectively a one-dimensional interaction, 
and all higher harmonics decouple from the impurity. 

In the presence of Rashba spin-orbit coupling, it is convenient to write 
the Kondo interaction in terms of the operators that create and annihilate 
Rashba quasi-particles. First we define the following operators, which are 
directly related to the m = and m = — 1 angular components of the Ak± 
operators respectively 




a fe ±T := y^J a d9A ± {k,6) (12) 
a k±l := d0e w A ± (k,6). (13) 

These satisfy the usual Fermionic anti-commutation relations {a\ bu , a-k'cu} = 
Sb c S^S(k— k'). By combining eq. (|10|) and the inversion of eq. one can 
show that 

tpkiio = -j= (ak+fj. + a k -fj,) ■ (14) 

Since ipk^o is the only electronic mode that couples to the impurity (cf. 
cq. (fTTjl ). we see that only the m = 0, — 1 modes of the Rashba quasi-particlc 
operators couple to the impurity and all other modes can be neglectec0. 



2 The fact that only the -1 mode couples and not the +1 mode can be traced back 
to the choice of global phase in the states that diagonalize the Rashba Hamiltonian Q. 
Of course, such a choice is arbitrary and does not effect the resulting analysis. 
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Hence, we may write the Kondo-Rashba Hamiltonian as a one dimensional 
model 



^ poo 

H = 79^2 / dk Yl £b ( fc ) a k afe ^ 
/>oo 

/ dkdk' \fk~k 1 V" S ■ a\ b a-^ak'cv (15) 
Jo , ' 



J V" 
4 (2tt) 3 



b,c=± 



This has the form of a two channel Kondo model with off-diagonal cou- 
plings, where each channel has a different dispersion relation and the cou- 
pling with the impurity is the same between each of the channels. In this 
expression, the Greek indices no longer refer to the spin of the electron but, 
instead, a combination of the spin and orbital angular degrees of freedom. 
This can be seen explicitly by writing the ak±^ in terms of the angular 
modes of the electrons (fTUj) 

Ofc±r = (^feto T iipkl+i) (16) 

ak±i = ^= (V'feio T iipkt-i) ■ (17) 

Since we are interested only in the low energy properties of the Kondo- 
Rashba model at this stage, we can focus only on the quasi-particle excita- 
tions close to the Fermi energy. Hence, we expand each of the dispersions 
to linear order about their respective Fermi momenta kF± 

e^fc) w e F + v F {k-k F ±). (18) 

We consider this approximation to be valid in some momenta range k £ 
(k F ± — A, k F ± + A) for each of the ± bands respectively, such that k R <C 
A -c k F . The expressions for the Fermi momenta and velocity as derived 
from eq. ([3]) are given by 



k F± = k° F Jl+^Tk R (19) 
£F 



v F 



d4 
dk 



1 + ^ (20) 
e F 



k F ± 

where we use a superscript to indicate values for the Fermi gas in the 
absence of spin-orbit coupling (i.e. k F = y/2me F , v F = ^2e F /m). An 
important point to note is that the Fermi velocities are the same for both 
the + and the — Rashba quasi-particles. 
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The kinetic energy of the Rashba quasi-particles (the first term in 
cq. ([15])) now reads 

1 f A 

Hr ~ ThrV VF ^ / dk ka lbn a kb» (21) 

where we have shifted the integration variable to be measured with respect 
to kp± in each of the two integrals in the summation respectively and 
simply relabelled the operators a(k+k F ±)±n l— > afe±^- The constant ef part 
of the linearized dispersion has also been dropped. 

We now similarly restrict the region of integration in the Kondo interac- 
tion (the second term in eq. (|15p ) and expand \fk and \fk! about either of 
kp± as appropriate. Since all higher order terms in k and k' are irrelevent 
we keep simply the constant term in each expansion. The result is that 
only the combination y/kp+ak+fj, + y/kp-ak-fi couples to the impurity. 
Hence, we define a new complete set of orthonormal operators 



V k F+ \kF- (v^U + V^U) (22) 
In terms of this new basis, our effective one dimensional model is 

H = J^f VF I A dkk { b l bk v + ~ b U^) 

+JCS J^T k ° F j\dkdk>S-bl^b k ,„. (24) 

If we neglect the free b quasi-particles that do not couple to the impurity, 
then this Hamiltonian is of the form of a single channel, single impurity 
low energy Kondo model with rescaled parameters Vp , given by eq. (|20[) . 
and J cff , given by 



As is necessary, this reduces to the original Kondo model with unsealed 
parameters in the limit sr — > 0. This is an example of the N channel model 
studied by Simon & Affleck [TS] where it was shown that such models can 
always be reduced to a one channel model provided that the Fermi velocities 
for all channels are the same (which is the case for the two Rashba bands 
as shown by eq. ([20])). 
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We can interpret the change in the Fermi velocity as a change in the 
effective mass of the quasi-particles that interact with the impurity. The 
corresponding change in the constant density of states for 2D electrons 
pa = m/2ir is then 

=^=. (26) 
1 + ff 

Multiplying ([25)) and (|26[) gives the result that 

J cS p = Jp . (27) 

Since the Kondo-Rashba model has been shown to reduce to the Kondo 
model with rescaled parameters, the Kondo temperature is the same as 
that of the ordinary Kondo model, eq. ©, with the above rescaled coupling 



T K = De '^ w ~p=De~~>. (28) 

That is, up to any possible changes to the energy cutoff D which, in this 
analysis, was implemented ad hoc, the behaviour of the Kondo-Rashba 
model is the same as that for the Kondo model and the value of the Kondo 
temperature is unchanged. 



3 Resistivity of the Two-Dimensional Model 

As we saw in the above section, any deviation in the low energy behaviour 
of the Kondo-Rashba model from that of the ordinary Kondo model (if 
there is any deviation at all) will be manifest as a change in the Kondo 
temperature via a change in the effective bandwidth cutoff D. In this 
section, we bring back the full two dimensional Kondo-Rashba model that 
includes the full energy range of the Rashba bands up to some high energy 
cutoff imposed by the finite range of the impurity interaction. In this way, 
band effects that were lost in the low energy linearization can be included in 
determining the change in the effective cutoff that appears in the expression 
for the Kondo temperature. Specifically, we calculate a physical quantity 
of the Kondo-Rashba model, namely, the resistivity, and compare it to the 
same quantity of the Kondo model, and interpret any difference as a change 
in this effective cutoff stemming from spin-orbit effects. The modification 
to the cutoff will then be incorporated into the expression for the Kondo 
temperature. 

As discussed in appendix [XJ a convenient choice for the finite range 
interaction with the localized Anderson impurity (in the absence of spin- 
orbit interaction) leads to a Kondo model of the form ([5]) with J kk i = 
Jh(k)h(k') where h(k) = 0(A — k), 9 being the Heaviside step function 
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and A = lix j a being the momentum scale associated with the range of the 
interaction a. 

We will calculate the conductivity using the full two dimensional Kondo- 
Rashba model H = Hr + Hk where the terms on the right hand side are 
given by @ and 0. Following the original calculation of Kondo pQ, we 
use the semi-classical Boltzmann theory of transport, appropriately gener- 
alized to the present case where the two bands are no longer degenerate, 
as described in appendix [Bl 

The resistivity is R(T) = l/a(T) where a is the conductivity given by 
equation ([50]) . specialized here to the case of the two dimensional ± Rashba 
bands 

'< r > = -y/(0„? ± «> 2 ^- < 2<J > 

Here, e is the electron charge, v*j: :— is the velocity of the a th band, f a 
is the Fermi-Dirac distribution for the a th band, and rf (k) is the inverse 
of the scattering rate of the a th band given by eq. (fMf 



T 



k 



(30) 

In this equation, Cj mp is the concentration of magnetic impurities (as- 
sumed to be dilute in order to neglect interactions between the impuri- 
ties themselves), 9' is the angle of k' as measured with respect to k, and 
Tka k'a' — (ka\T\k'a'} is the T-matrix, i.e. matrix elements of the operator 

T = H K + H K G+H K + H K G+H K G+H K + ■■■ (31) 

written in the basis of eigenstates \ka) of Hr with = l/(e — Hr + iif) 
being the retarded Greens function. The angle brackets indicate a free spin 
averaging over the impurity operators in the T-matrix. 

If we consider J <C £_f then we can interpret eq. pip as a perturbation 
expansion and compute it to 0(J 2 ) using the full Kondo-Rashba interac- 
tion , as denoted diagramatically in Fig. [TJ Keeping terms to third order 
upon squaring and spin averaging these matrix elements results in 

Tu \ I ) = ( \Ti 



fe+,fe'+| / — yl-tfe-.fe'- 

J \h(k)h(k')f (l + sin 2 



2 



l--(g++g +C.C.) 



(32) 
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q+ 




+ 




Figure 1: The diagrammatic expansion of the second order contributions 
to the T matrix element Ty + k + for the finite range Kondo-Rashba model. 
The solid line represents the propagator for a Rashba quasi-particle, the 
dashed line represents the impurity, and time flows from left to right. The 
internal momentum q is summed over. After interacting with one of the 
Rashba quasi-particles the impurity is no longer in a spin eigenstate and 
so the spin states of the impurity are not labelled here. 



and 



T, 



fe-,fe'+ 



'-k+,k'- 



J l\h{k)h{k'f (l + cos 2 I) 



l--(g+ + g +C.C.) 



(33) 



In these equations, 0' is the angle between the two momentum vectors k 
and fe', c.c. represents the complex conjugate of the preceding terms, and 
we have defined 



V 



(27T)2 



d 2 q\h(q)\ 



2 / ( £ g = ) 2 

- ef +irj 



(34) 



which comes from summing over the intermediate momentum state and is 
non-divergent due to the momentum cut-off coming from the finite size of 
the interaction via the function h(q). 

In order to continue with this calculation, we make some assumptions 
regarding the various energy scales that are present in this system. Let 
D := A 2 /2m be the energy scale associated with the high-energy cutoff 
that arises from the finite range of the interaction and assume that D > ep. 
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Relative to the Fermi energy, the lower cutoff is approximately £p (the 
distance between the Fermi energy and the bottom of the band) and the 
upper cutoff is approximately D — Ep- We have already assumed that 
J <C £f by truncating the above perturbation expansion. Since the spin- 
orbit coupling is relatively weak in all physical systems that are studied, 
we also assume that er <C £f and er <C D. We are interested in the 
low temperature behaviour and can therefore assume that T <C £f and 
T <C {D — ep)- This last assumption allows us to approximate derivatives 
of the Fermi-Dirac distribution 

%~-5{e-ef) (35) 
as 

except in terms that diverge as T — > 0. Such an approximation, as was used 
in Kondo's original calculation [T], is valid since the dominant behaviour 
comes from the diverging logarithm, as will be seen. 

Substituting eqs. ([32]) and ([33)) into (|30[) results in the scattering rates 



rt{k)- 8 M 2 



^[l--(g + +g-+c.c.) . (36) 



The function takes different forms depending on which range of k space 
one is in with respect to kR and A. However, we are only concerned with 
the value of this function in the vicinity of the Fermi wave vectors since 
Ty is eventually multiplied by a derivative of the Fermi-Dirac distribution 
in eq. (|29|) . Hence, we can take 



which is strictly valid for 2kn < k < A — 2fc^. 

To evaluate and the integral in the expression for the conductivity, 
it is convenient to convert them to integrals over energy using the density 
of states which, for the Rashba bands, are 



1- ' £R 



p + (e) = Po e(e) 

V £ + £ R 

P~{e) = po&{£ + £r) 



£ + £r \ V £ + £ R 



(38) 
(39) 



where, as above, po = m/2ir is the two dimensional density of states for a 
Fermi gas. Under the approximations discussed above, inverting eq. (|36p 
to leading order in J and £r and substituting into eq. (|29|) yields 

4e 2 (e F + §£fl) / J [ D , / + , - , x df S 



" T »° TO '-'l (40) 
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The difference in the iimits of integration between the two bands has been 
ignored on account of the presence of df/de. 

However, the integrals in g are evaluated by integrating by parts. The 
integrands are evaluated at their respective limits and so the differences 
between the two bands will become manifest. Writing 
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de'p ± {e') /i(fc ± (e')) 



e — e 



17] 



(41) 



where fc ± (e / ) is the inversion of , integrating by parts, substituting into 
cq. (|40p . and inverting to leading order in J yields an expression for the 
resistivity 



R(T) 



a(T) 4e 2 (e F 



37r 2 c imp J 2 po 



J 



(I + PoK) 



The diverging logarithm comes from the term 



a,b=± 



de 



cfe'ln 



e F (e b A - £f) 



dj_d[_ 

de de' 



(42) 



(43) 



and all of the terms arising from the integral over the square root in the 
density of states have been collected in the term 



K 



SF+6R 




arctanh-i 



£R 



(44) 



Again, we have used the assumptions and approximations discussed above 
except in the case of the log term which diverges as T — > 0. 

We can extract the temperature dependence by making the change of 



variables 
expression 



r(0 



£f)/T in the above integral to obtain the final 



R(T) 



37r 2 c im pJ 2 /?o 
4e 2 (e F + §£«) 



-po In 



4( £ A - £ f){£ a ~£f) 



K + C 



(45) 

where C is simply a finite number. As en — > 0, this reduces to the classic 
resistivity of the Kondo model 



R°(T) = 37r2 Cim P j2po 



4e 2 £F 



1 • 1 ln 77T ' T 



(46) 



with D° := ^Je F {D — e F ) being the UV cutoff in the absence of spin-orbit 
coupling. 
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3.1 Interpretation of spin-orbit effects 



The analysis of the effective one dimensional model showed us that any 
change in the low-energy behaviour of the Kondo-Rashba system should 
be interpretted as a change in the cutoff. Since we are considering the case 
where er<^£f and er -C D, we can simplify the above expression for the 
resistivity eq. (|45|) . by keeping only terms up to first order in sr/ef and 
sr/D 



R(T) 



-'imp 



4e 2 (e f + \e r 



1 -- 2J ^ ln £ + T 



(47) 



where we have interpretted all of the changes as giving rise to an effective 
cutoff 

Given that we had no way of determining the change in the cutoff 
from the low energy one dimensional model, we propose that the change 
in the cutoff determined from the calculation of the resistivity may give 
us a further clue since it incorporates the full, non-linearized bands and a 
physical "real" cutoff. In the presence of Rashba spin-orbit coupling, the 
Kondo temperature would be 



^ ( e r D{D-2ef) 

~ "*° - T * ( 1 + T F (D - sfY 



(49) 



where Tg- is the Kondo temperature using the cutoff D°. 

In the case where the magnetic impurity is extremely localized, D ^> ef, 
this expression reduces to 

T K « T° + (50) 

which is independent of the details of the cutoff. The generality and utility 
of these expressions will be discussed in the following section. 



4 Discussion 

The analysis presented above shows that the low energy behaviour of the 
Kondo model, that is, the formation of a singlet that screens the magnetic 
impurity, persists in the presence of a Rashba spin orbit coupling even 
though such an interaction breaks the SU (2) symmetry so that spin states 
are no longer eigenstates of the Hamiltonian. The effective one dimensional 
model for the Kondo-Rashba system has the form of a single impurity, single 
channel Kondo model with an effective mass and exchange interaction that 
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is rescaled as a function of the strength of the spin-orbit coupling, eqs. ([^5]) 
and ([25]) . However, this rescaling combines in such away as to be cancelled 
out in the expression for the Kondo temperature, which is the only energy 
scale on which the low energy behaviour of this system depends. However, 
to derive such a low energy model, one must impose an ad hoc cutoff 
which may be modified in the presence of spin-orbit coupling, though it is 
impossible to derive such a change of the effective cutoff in this low energy 
formalism. 

To get a better sense of how the cutoff may be modified, the resistivity 
of the full two dimensional Kondo-Rashba model, defined over the whole 
energy range up to a high energy cutoff, D, arising from the finite size of the 
interaction, was calculated and changes wrought by the Rashba coupling 
interpretted as giving rise to an effective cutoff, eq. (f48|) . Using such a cutoff 
in the expression for the Kondo temperature leads to the prediction that 
the presence of Rashba spin-orbit coupling will have the effect of slightly 
increasing or decreasing (depending on whether D is greater than or less 
than 2ef respectively) the Kondo temperature as a linear function oien/eF 
for small values of this dimensionless parameter. 

In the calculation of the resistivity, all of the modifications that arose 
due to the spin-orbit coupling came from the changes in the band structure 
(as was manifest in the densities of state /5 ± (e)) and not from the modifi- 
cations to the Kondo coupling, despite the apparent complicated nature of 
the interaction (J7J). Hence, it is not surprising that the low energy effective 
model, which required us to linearize the energy bands about the Fermi 
energy and so neglect the changes in the band structure (since both bands 
have the same Fermi velocity), failed to show any change in the physical 
behaviour of the system. 

Yet, it was still necessary to impose a cutoff in the full two dimensional 
model by introducing a finite range to the hybridization between the impu- 
rity and the conduction electrons, the precise form of which (characterized 
by the envelope function h(k)) was chosen for simplicity. Hence, any results 
that depend on the nature of the hybridization may be of questionable use. 
In particular, the general expression for the Kondo temperature, eq. (|49[) . 
depends on the energy cutoff D. However, in the case of a highly localized 
impurity, D ^> sf and the universal result stated in eq. (|50|) is obtained. It 
is in such a localized regime that these results are expected to be of greatest 
utility since they are independent of the details of the cutoff mechanism0 

3 It should be noted that the universal expression for the Kondo temperature l|50j l 
still assumes that the hybridization between the local spin and the conduction electrons 
is circularly symmetric. 
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4.1 Application to Experiment 

The simple models used to derive these results were chosen so that the 
effects of the spin-orbit coupling would be as transparent as possible. One 
may now inquire as to their applicability to real two dimensional systems 
that are of experimental interest. In particular, we consider two broad pos- 
sibilities: that of quantum dots constructed in semi-conductor heterostruc- 
tures and that of magnetic atoms placed on a clean metallic surface. 

In the case of quantum dots, we conclude that the above analysis cannot 
apply for two reasons. For one, the imposition of rotational invariance in 
the plain of the electron gas is clearly not present in the complicated gated 
structures that are manufactured in the laboratory. In our calculation of 
the resistivity, rotational invariance played a crucial role in maintaining 
the Kondo effect (as manifest by the logarithm term) in the presence of 
spin-orbit interactions. To be more precise, the rotational invariance is 
the reason why the integrations over an intermediate -I- and — Rashba 
quasi-particle in the perturbative calculation of the T-matrix (as denoted 
diagrammatically in Fig. [T] and represented analytically by the terms g + 
and g~ respectively in eqs. (|32"j) and {33])) occur with equal weighting in the 
final T-matrix expression; the difference between the interactions involving 
+ and — particles, as shown in the Hamiltonian ([7]), all occur with a factor 
of e ±l0 in the T-matrix and so integrate to zero only when you integrate 
over all angles from to 2ir. This is reminiscent of the results of [21] 
wherein a magnetic impurity embedded in a host metal with spin-orbit 
interaction was studied. It was shown that the Kondo effect is supressed 
only when the impurity is close to a boundary and, hence, the rotational 
isotropy of the system is broken. 

The second reason why our analysis fails in the case of semi-conductor 
quantum dots is the relative strength of the energy scales. In particular, if 
one assumes that the distance scale of the finite range, a, is of order the 
radius of a typical quantum dot (say a ~ 50 nm [3]), then this leads to an en- 
ergy cutoff D = A 2 /2m = 2n 2 / (ma 2 ) w 0.008 eV. This is generally smaller 
than the typical values of the Fermi energy in the two dimensional electron 
gas of the semi-conductor heterostructure (typically Ef ~ 0.03 eV [18]), 
completely negating our necessary assumption that D > ep- This suggests 
that a more accurate model is required to study the Kondo effect with 
Rashba coupling in quantum dots, some examples of which have already 
been studied [H EH HI EDS • 

It is more likely that the results obtained in this paper are applicable 
to systems of magnetic atoms placed on the surface of a metal. There, the 
rotational isotropy is manifest (assuming that the surface is "clean" ) , and 
the dispersion relation of the surface conduction electrons more closely ap- 
proximated by the two Rashba bands ([3]) (as shown in various experiments, 



1G 



e.g. those described in [10). Furthermore, the Fermi energy is generally 
much larger for the surface state electrons than in the semiconductor two 
dimensional electron gas (ep ~ 0.5 eV for typical metallic surfaces) mak- 
ing it much more likely that the length scale over which the magnetic atom 
hybridizes with the conduction electrons is such that D > Ep. Indeed, 
assuming that such a hybridization takes place over the observed size of 
the impurity atom, a radius of a ss 7 A [10] . then D w lleV > Ep. 

As a result, one could imagine detecting a change in the Kondo temper- 
ature as a result of the Rashba spin-orbit coupling by measuring the Kondo 
temperature of the same magnetic atom on different metallic surfaces with 
different strengths of the Rashba coupling. This is likely to be terribly 
difficult to detect given that the change due to spin-orbit coupling is typ- 
ically very small, being linear in er/ep (eji/ep » 0.003 for the case of an 
Au(lll) surface and immeasurable in Cu(lll) and Ag(lll) [lQl US])- Such 
an effect is likely to be overwhelmed by other differences between the two 
systems. A clearer signature of this effect could be seen if one were able to 
tune the Rashba coupling independently from the Fermi energy, something 
that has been demonstrated in semi-conductor heterostructures [22] but 
not for the coupling in metallic surface electrons. Hence, for the purposes 
of observing this effect, it may be worth pursuing a more realistic model of 
a quantum dot system in such a two dimensional electron gas where one 
has the ability to tune the Rashba coupling. 

As was demonstrated in the presence of time reversal symmetry 
plays a fundamental role in the survival of the Kondo effect in the presence 
of spin-orbit interactions. Another interesting question to pursue is how 
the presence of a magnetic field, which breaks time reversal invariance, can 
be incorporated into the models introduced above and how such a field 
modifies the results. 
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A Derivation of the Kondo model from a fi- 
nite range Anderson impurity model 

Here, we review the Schrieffer- Wolff transformation of an Anderson impu- 
rity model without spin-orbit coupling. The purpose is to highlight how a fi- 
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nite range hybridization between the impurity state and the non-interacting 
conduction electrons manifests itself in the resulting Kondo interaction. It 
is that Kondo interaction that is then used as input into the two dimen- 
sional Kondo-Rashba model studied in this paper. 

The Kondo model used in this paper can be thought of as deriving 
from the Anderson impurity model [T7] of an impurity electronic state 
annihilated by interacting with a sea of conduction electrons annihilated 
by ^kfi- The Hamiltonian for this model can be written as 

H = j-^ J d 2 ke k y{^ kt , + eadldf, + Un dr n dl + H hyh (51) 

where n djl :— dj^d^y The hybridization between the impurity and conduc- 
tion electrons is chosen to take the following form 

#hyb = J d 2 k (h(k)dl* kfi + h.c.) . (52) 

Here, V is the two dimensional volume of the system, h.c. denotes the Her- 
mitian conjugate of the preceeding term and h(k) is the Fourier transform 
of an "envelope function" h(x) that describes the extent over which the 
impurity interacts with the conduction electrons. This envelope function 
has the property that its maximum is at the origin, that it only depends 
on the distance from the origin (i.e. it is circularly symmetric), and that it 
decays to zero after some characteristic length scale a. 

One can perform a canonical transformation, the Schrieffer- Wolff trans- 
formation [20], perturbatively in t/U, to obtain the s-d or Kondo model 

H= hJ rf2fe **l^ + / d 2 kd 2 k>J kk ,S ■ *l^V h , v (53) 
where S = dj^^^-d^ is the operator of the impurity electron and 
J kk , =t 2 h*(k)h(k') 
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U + e d - £fe' e k - e d 



(54) 



A potential scattering term is also generated but is neglected in our analysis 
as it does not play a crucial role in determining spin-orbit effects. If we 
are deep in the so-called local-moment regime where the singly and doubly 
occupied impurity levels are equally spaced below and above the Fermi 
energy {i.e. e k = e k i ~ £f and e d ~ £f — U/2) then we can simplify the 
above to 

Jkk' w —h*(k)h(k') ee Jh*(k)h(k'). (55) 
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The traditional choice in most theoretical analysis is to take h(k) = 1, 
which is the Fourier transform of a delta function in position space and 
corresponds to having a point-like interaction. This is the form of the 
envelope we will use in § 12.21 However, in § [31 it will be necessary to make 
the range of interaction finite in order to provide a more physical ultra- 
violet cutoff. So, for ease of computation, we take h(k) = 6(A — k) which 
corresponds to a position space envelope function 

h(x) = -MAx) (56) 
x 

with J v {z) being the Bessel functions of the 1st kind. Such a function has 
an effective range of a = 2tt/A. 

B Boltzmann theory of transport for non- 
degenerate bands 

In this appendix, we generalize the standard Boltzmann theory of trans- 
port, as described in [2j, for example, to the case where there are multiple 
bands, eJJ, that are not necessarily degenerate. Consider such a system in 
the presence of a weak static electric field E. This will give rise to a distri- 
bution function, fg(k), for each band. The time rate of change is assumed 
to satisfy the relaxation hypothesis 

df a E (k) /g(fc) - f a (k) 

dt rf(fc) [ ' 

where f a (k) = /(£?), the Fermi-Dirac distribution. Assuming that such 
time rate of change only comes about via the change in momentum k as 
a function of time, then we can equate the above to (dk/dt) ■ V/Jj(fe) = 
— eE ■ V/^(fc). Working to first order in E yields 

f a E (k) « f a (k) + eT -(k)^^-vt ■ E (58) 

with vf. := Ve£ being the band velocity. 

Substituting this into the expression for the current 

J = (t(T)E = -ej ' ^ £ f a E {k)vt (59) 
and averaging over all directions of space gives us 
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To calculate the relaxation time, we introduce the scattering rate be- 
tween the state \ka) and \k a') given by Fermi's golden rule 

W ka ^ a , ^2n(\ [ T ka ^ a ,\ 2 )8{e a k ^ei,) (61) 

where T kak > a i is the T-matrix as introduced in § [3] 

The rate of change of the distribution function f%(k) due to impurity 
interactions can then be obtained by summing over all possible transitions 
of the state \ka): 



df a E (k) 
dt 



-•imp 



d d k' 



E 



-W k - a , M fi{k'){l- f%{k)) 



= -2nc imp f ^- d Y^(\T ka , k , a ,\ 2 ) (f%{k) - f a E (k')) 6(e a k - 4), 

a' 

Cimp being the concentration of impurities. Equating this to — eE ■ V/|(fc) 
and using eq. (|58p results in the expression 



E = 2irc imp T?(k) 



d d k' 



(2tt) 



d Z^i 



T kaMa ,\ 2 )(vl-vi^-E8(e a k 



ez, 



(62) 

Specializing now to two dimensions, we assume that vf, is at most linear 
in k (which is true for the Rashba bands studied in this paper), that the 
T-matrix depends only on the angle between its two input momentum, and 
that the T-matrix is symmetric in that angle. These allow us to make the 
following substitution within the integrand above 



vi, ■ E ->vt,Ecos0cos0' 



V-vl-EcosO' 



(63) 



where 9 is the angle between k and E and 9' is the angle between k and 
k' . The reason this works is that the sine term that would also be included 
in the expression for the dot product vanishes after integrating over the 
angular variable. Using this substitution we can extract an expression for 
the scattering rate 



d d k' 

{2ir) d 



E (l^^'l 2 ) ( 1 - ^ cos9 ' ) § ( s t ~ 4) 



(64) 



This can then be used in eq. (|60j) to evaluate the conductivity. 
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